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Abstract 

We investigate the convergence properties of the cluster expansion of equal- 
time Green functions in scalar theories with quartic self-coupling in (0 + 1), 
(1 + 1), and (2 + 1) space-time dimensions. The computations are carried out 
within the equal-time correlation dynamics approach, which consists in a closed 
set of coupled equations of motion for connected Green functions as obtained 
by a truncation of the BBGKY hierarchy. We find that the cluster expansion 
shows good convergence as long as the system is in a localized state (single phase 
configuration) and that it breaks down in a non-localized state (two phase con- 
figuration), as one would naively expect. Furthermore, in the case of dynamical 
calculations with a time dependent Hamiltonian for the evaluation of the effective 
potential we find two timescales determining the adiabaticity of the propagation; 
these are the time required for adiabaticity in the single phase region and the 
time required for tunneling into the non-localized lowest energy state in the two 
phase region. Our calculations show a good convergence for the effective poten- 
tials in (1 + 1) and (2 + 1) space-time dimensions since tunneling is suppressed 
in higher space-time dimensions. 



PACS: 11.30.Qc; 11.90.+t. 



* supported by DFG, Forschungszentrum Jiilich and GSI Darmstadt 



1 



1 Introduction 



The many-body problem of quantum field theory, in spite of great efforts in the past, 
still remains unsolved and there is a need for genuine nonperturbative methods. In 
M we have proposed a connected (equal-time) Green function approach for SU(N) 
gauge theories which might provide valuable insight into the low energy QCD problem. 
However, the possible truncation schemes in the order of the connected Green functions 
will be limited for practical purposes and one thus needs intrinsic criteria that allow 
for a judgement of the convergence properties of the approach especially in those cases, 
where exact solutions are not available. 

Before adressing the SU(N) Yang-Mills problem, we analyze the convergence prop- 
erties of the correlation dynamical approach for the scalar quantum field theory with 
A$ 4 self-interaction, which has become an important theoretical laboratory for test- 
ing non-perturbative methods. This is partly due to the fact that $ 4 -theory provides 
an easy framework for studying different scenarios in nonperturbative renormalization 
when going from (0 + 1) to (1 + 1) §-0, (2 + 1) 0-Q and (3 + 1) @-0 space-time 
dimensions. 

In two preceding articles |15[ [l(| we have investigated ground state symmetry break- 
ing and the effective potential in <3> 4 -theory in (1 + 1) and (2 + 1) space-time dimensions 
using the equal-time correlation dynamics method (in (o?+l)-dimensional $ 4 -theory de- 
noted as &j +1 CD), which consists in a coupled set of equations of motion for connected 
equal-time Green functions. This coupled set of equations is obtained by inserting the 
cluster expansion, i.e. the expansion of full Green functions in terms of connected 
Green functions, into the equations of the BBGKY hierarchy and by neglecting all 
connected n-Point functions with n > N for some given N. 

The derivation of the correlation dynamical equations for $ 4 -theory has been pre- 
sented in detail in |[|, [16]; the resulting equations in (1 + 1) dimensions are given by 
Eqs. (20)-(31) in appendix A of [IK] and in (2 + 1) dimensions by Eqs. (19)-(21) in 



Sect. 2 as well as Eqs. (28)-(36) in the appendix of JT^]. The equations in (0 + 1) 
dimensions can easily be obtained from those in higher dimensions by discarding the 
dependence on the spatial coordinates and by omitting all mass counterterms. For 
a more general discussion of correlation dynamics in case of nonrelativistic fermionic 
systems we refer the reader to refs. |T7|, [18], [HJ f20] . 

The calculations in JIB], [0| aimed at the evaluation of the effective potential V e ff($o), 
which is given by the minimum of the energy expectation value in the subspace of states 
with a fixed expectation value of the scalar field ($) = $ - Since the CD-equations 
describe the (equal-time) evolution of the system for a given initial condition, a direct 
(variational) evaluation of the effect iv potential within this approach would amount 
to minimizing the energy in the manifold of stationary solutions with the constraint 
(<&) = $o- In practice, however, this turns out to be impossible because the energy 
is not bounded from below without additional constraints, since the solution manifold 
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contains unphysical Green functions. Although e.g. for the (0 + l)-dimensional case in 
the 2-point approximation a sufficient constraint can be derived from the uncertainty 
relation, it is not yet clear what constraints for higher order approximations have to 
be imposed. 

The computation of the effective potential and the corresponding lowest energy state 
Green functions within our approach is therefore carried out by means of a dynamical 
calculation using the Gell-Mann and Low theorem, i.e. by adiabatically changing the 
parameters of the theory in time from a configuration with a well known lowest energy 



state solution to the actual configuration of interest as described in [To, 16). Whereas 



this many-body scheme was found to converge quite rapidly in case of nonrelativistic 



nuclear physics problems [20|, it is not clear if this also holds for the quantum field 
theoretical case. We thus have to explore the convergence of the cluster expansion and 
its dependence on the parameters of the theory, i.e. the coupling A and the expectation 
value of the scalar field (<&) = $o- 

Before investigating $ 4 -theory in (1 + 1) and (2 + 1) space-time dimensions, we 
turn to the simpler case of (0 + 1) space-time dimensions, which is equivalent to the 
quantum mechanical anharmonic oscillator in case of a positive squared mass and to 
the quantum mechanical double-well potential in case of a negative squared mass. Since 
in (0 + 1) dimensions we have access to the exact lowest energy state solution for given 
(x) = x , we here can also directly control the validity of the correlation dynamical 
solution. 

This article is organized as follows: In Sect. we investigate the (0 + l)-dimensional 
case in full detail with respect to the problem of degenerate vacua and quantum tun- 
neling, whereas in Sect. |3| we present our numerical results for the (1 + 1)- and the 
(2 + l)-dimensional case. Sect. [| concludes the present study with a brief summary 
of our results, while technical details of the cluster expansion as well as the GEP 
(Gaussian effective potential) approximation are given in appendices [A] to |C]. 



2 Convergence of the cluster expansion in (0 + 1) 
dimensions 

In (0 + 1) space-time dimensions and with the identification $(£) — > x(t), TL(t) — > p(t), 
the Lagrangian of $ 4 -theory with positive squared mass m 2 (anharmonic oscillator) 
reads 

Li = ~ - \\^ , (1) 

while the Hamiltonian is given by 

Hl = i p 2 + I m 2 x 2 + i Ax 4 (2) 

T 2 4 v ' 
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with the usual equal-time commutation relation [x(i),p(i)] = i. 

In case of a negative squared mass m 2 = —^i 2 (double-well potential) we write the 
Lagrangian and the Hamiltonian as 

L2 = kd t *f + - V - £ = W - iA(x 2 - (3) 



4A 2 V ' 4 v A 



and 



1 \ u 2 M 



2 



^2 = -P 2 + ~A(x 2 - ^-) 2 . (4) 

Since we are aiming at the effective potential, i.e. the properties of the system not 
only as a function of the coupling A, but also as a function of (x) = x , we explicitly 
split off the vacuum expectation value by writing 

x = x + x , (x) = ; p = p (5) 

and treat xo as a constant classical background field which can be varied as an external 
parameter. The explicit $q +1 CZ) equations then are obtained in the usual way, i.e. by 
inserting the cluster expansion (cf. appendix [S]) into the coupled set of equations of 
motion for the n-point Green functions (xx), (pa;), (pp), (xxx), (pxx), (ppx), (ppp), 
etc. The resulting set of coupled equations for the connected Green functions (xx) c , 
(px) c , (pp)c, then can be closed by neglecting connected Green functions above some 
order n > N, where N is an integer > 2. Since the resulting set of equations is quite 



lengthy and in case of N = 4 has been presented in refs. Hl5 , [T6|1 for (1 + 1) and (2 + 1) 
dimensions, we omit an explicit representation here. 

In order to investigate the convergence of the cluster expansion, we compute the 
quantities 

l(x n ) c /(x")|, (6) 

where (-) c stands for the connected part of the expectation value (•) (cf. appendix 
for n = 2, 3, 4, 5, 6 in the exact lowest energy state solution as well as in the 
^o+i^^ approach using a 6-point truncation scheme, i.e. including the connected n- 
point functions with n < 6. In line with the notations in [13. W\ this set of equations 
is denoted by §* +1 CD(2, 3, 4, 5, 6). 

In our analysis we consider two quite distinct cases; i.e. an extreme double- 
well potential with A/(4/i 3 ) = 0.0158 and an anharmonic oscillator potential with 
A/(4m 3 ) = 10. The values (^) of the exact solution have been obtained by diago- 
nalizing the Hamiltonian of the system with an additional external source (cf. @) 
in a sufficiently large set of basis states. The §q +1 CD(2, 3, 4, 5, 6)-results have been 
obtained in the usual way by exploiting the Gell-Mann and Low theorem. We start 
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with the GEP solution |7|, [l^, [21] (cf. appendix |C|) and switch on the residual inter- 
action terms adiabatically. In case of the anharmonic oscillator this corresponds to a 
time-dependent Hamiltonian of the form 



H (t) = y + + 




x 2 + -Ax' 
4 



(7) 



where the residual coupling g(t) runs from zero to A, M is the effective mass of the GEP- 
solution and x is fixed. We recall that in the correlation dynamical calculation the 
1-point functions (x) and (p) are set equal to zero according to (JD, and for each value 
of x an individual time dependent calculation has to be carried out. We have used 
a linear time dependence of the form g(t)/4fi 3 = f3t or g(t)/4m 3 = (3t, respectively, 
where (3 in each case has been chosen sufficiently small such that the system shows 
adiabatic convergence as long as xo is outside the domain where the exact solution 
indicates a breakdown of the cluster expansion. Our numerical results for (§) in the 
lowest energy state with given expectation value x are shown in Fig. |as a function 
of x for both potentials. The explicit results from the exact solution are displayed on 
the l.h.s. whereas the r.h.s. shows the corresponding quantities in $q +1 CD(2, 3, 4, 5, 6) 
approximation. We note that the sharp peaks for certain values of x /y/m or x /-v/p 
in Fig. |] result from a change of sign of (x n ) c / (x n ). 

In the case of the double-well potential the exact solution (upper left part of Fig. 
[1]) yields a dramatic sudden change in the relative importance of the connected Green 
functions at a value of Xq" 4 / \fji ~ 3.8. Above x c ^ %t j ^fji the relative importance of the 
connected Green functions decreases by at least one order of magnitude as one goes 
from nton + 1, implying an excellent convergence of the cluster expansion. In contrast 
to that, in the region below x^ lt / ^JJi the full Green functions are dominated by their 
connected parts, i.e. the system is dominated by fluctuations. 

The corresponding result in $q +1 C-D(2, 3, 4, 5, 6)-approximation (upper right part 
of Fig. [I]) above x^ lt / ^ffi shows nearly the same result as the exact solution. This 
has been expected due to the excellent convergence of the cluster expansion in this 
region. The sudden rise in the relative importance of the higher order connected Green 
functions in going from x / \J~ji > Xq"* / JJi to x / \fjl < XQ r **/^/Ju, however, is smeared 
out over an extended region for the finite value of (3 taken in the calculation. However, 
it becomes evident that as soon as the connected Green functions dominate the full 
Green functions (xq/ \JJi ~ 2.84), the method of time dependently switching on the 
residual interaction within the $q +1 CZ)(2, 3, 4, 5, 6)-approximation breaks down. 

In contrast to the results for the double-well potential, the ratios (§) for the anhar- 
monic oscillator (lower two pictures in Fig. [I]) demonstrate a good agreement between 
the exact solution and the $q +1 CD(2, 3, 4, 5, 6)-solution for all Xq/^//!. The relative 
importance of the connected Green functions increases with decreasing x /y/JI, but 
even at xq/-JJI = the higher order full Green functions can in a good approxima- 
tion be expressed by their disconnected parts (note that (x 2 ) c / (x 2 ) — > 1 as x — > by 
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definition) . 

In order to understand the numerical results for the double well potential in Fig. [I] 
in a more qualitative way, we now examine the wave functions obtained by the exact 
solution for the system. We first note that instead of finding the lowest energy state in 
the subspace of all states with given (x) = x , one can equivalently introduce a source 
term into the Hamiltonian by writing 

Hi = H 2 - (jtJx (8) 

and then find the ground state of the system with respect to the whole Hilbert space. 
The resulting state, which has some vacuum expectation value x ( J) depending on the 
source J, is also the lowest energy state in the subspace with given x (J). 

In Fig. ^| we display the absolute square of the ground state wavefunction ^(x) and 
the wavefunction of the first excited state x I / i(x) of the Hamiltonian (|]) with A/(4^i 3 ) = 
0.0158 as a function of f{J) = xo(J)/y / /I and x/yZ/I where both wavefunctions have 
been obtained via diagonalization of (§). At f{J) = 0, |\l/o(x) 2 | and |\l/i(x) 2 | have 
approximately the same shape with maxima at ±x mm . We note that in the limit 
A/ (4/x 3 ) — > the double- well potential develops an infinitely high barrier at x = 0. This 
leads to a spectrum of the Hamiltonian for J = which consists of pairs of degenerate 
Eigenvalues. Due to the infinite barrier at x = there is no more tunneling between the 
2 orthogonal states spanning the Eigenspace corresponding to each Eigenvalue. The 
lowest energy states of the system then are of the form 

|#o> =a\L ) + b\R ) (9) 

with \L ) and \Rq) denoting the mutually orthogonal lowest energy solutions with 
(L |x|L ) = -x min and (R \x\R ) = +x min , i.e. the optimally localized "left" and 
"right" states which can be constructed from an appropriate superposition of the de- 
generate ground states. We then obtain 

(#oW*o) = x™ n (|6| 2 - \a\ 2 ) e [_ x mi ",x™ n ] (10) 

with (Lo|x|i?o) = 0; states of the form (|9]) thus restore the convex shape of the effective 
potential by making it flat between — x mm and x min . As soon as J ^ 0, the system 
will energetically prefer to shift its ground state wavefunction completely to one side of 
the potential barrier, depending on the sign of J. In other words, the system with an 
infinitely high potential barrier undergoes a phase transition from the state |Lo) to the 
state \Ro) as J passes through zero. Here we use the term "phase transition" for any 
discontinuous behaviour of the lowest energy solution with respect to the parameters 
of the Hamiltonian. 

If A/(4/z 3 ) assumes a finite (but still small) value as in case of Fig. 0, the pairwise 
degeneracy of Eigenvalues disappears due to tunneling - which is no longer forbidden 
for finite A - but the spectrum still consists of pairs with a small energy gap. The 
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separation of the lowest Eigenvalues AE = E\ — Eq can be taken as a measure of the 
inverse tunneling time. The region between — x mm and x mm then can be covered by a 
very small variation of J, i.e. the above obtained "phase transition" from \Lq) to \Rq) at 
J = in case of a finite coupling is smeared out into a sharp crossover. As can be seen 
from Fig. |2], the ground state wavefunctions are (to a very good approximation) still 
of the form (|9]), with only slight modifications due to tunneling. The first excited state 
of the system then is given by the Eigenstate corresponding to the higher Eigenvalue 
in the lowest lying pair, i.e. approximately by 

\* 1 ) = b*\L )-a*\Il o ). (11) 

In Fig. [|, this can be seen from the fact that the relative strengths of the left and right 
maxima of ^(x) 2 ! and |\l/o( x ) 2 | depend on f(J) in exactly the opposite way. 

For |/(J)| > x mm / ' JJL the system suddenly changes its behaviour and the ground 
state and the first excited state are localized states, which to a good approximation can 
be described as the ground state and the first excited state of the harmonic oscillator 
potential obtained in second order Taylor expansion of the classical potential around 
its minima, as can also be observed in Fig. |2]. This corresponds to the single phase 
localized configuration at finite J in the limiting case of the infinitely high potential 
barrier. 

We arrive at the conclusion that for a finite, but small value of A/(4/x 3 ) - although 
the system does not exhibit true ground state symmetry breaking due to tunneling 
-, it will still qualitatively behave like a system which is in a (non-localized) two 
phase configuration for |xo| < x mm and in a (localized) single phase configuration for 
|x | > x min . 

It is now straightforward to show analytically, that the cluster expansion has to 
break down in a two phase configuration according to (0) as soon as it shows good 
convergence in both optimally localized single phase configurations. Furthermore, in 
a field theory with an arbitrary space-time dimensionality, the cluster decomposition 
property of the connected Green functions, 

(TflXTl \ 
i=l,-,p \ X i ~ Vj\ -* 00 > (I 2 ) 
j=p+l,-,n ) 

is lost in a state constructed according to (||) whenever \L ) and \R ) are the optimally 
localized states in case of a system with ground state symmetry breaking. For the (0+1) 
dimensional system this has no direct implication for equal-time Green functions since 
it only applies to the time arguments, e.g. in a time ordered Green function. 

Taking into account the above considerations, we can now identify x mm , the expec- 
tation value of x in the optimally localized right state, with Xq™*, for which (0) (cf. the 
upper left part of Fig. [1]) shows a relatively sharp rise. In the correlation dynamical 
calculation (displayed in the upper right part of Fig. IT]), the sharp rise is smeared out. 
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We also observe that the region, in which the system is dominated by fluctuations, 
only extends up to x / ^fji ~ 3.3 instead of xo/y/JJ = Xq"* « 3.8 for the exact solution. 

In order to get an explanation for this behaviour, we show in Fig. ^ a qualitative 
picture of the external source J as a function of xo for 3 different cases. In case (1) 
(dashed line), there is no ground state symmetry breaking and for each value of J there 
is a unique value xo(J). In case (2), furthermore, we have the situation corresponding 
to the convex shape of the effective potential given by the construction of nonlocal 
states (H) in case of ground state symmetry breaking, where for any nonzero value of 
J there is a unique value x (J) and for J = all values between — x mm and x mm are 
accessible. This corresponds to the field theoretical Maxwell construction [p2]| , where 
the two phases, given by the left and the right optimally localized state in (||) " coexist" . 
In case (3) we still have ground state symmetry breaking, but the curve for J(x ) is 
analytically continued into the region, where the field theoretical Maxwell construction 
takes place, implying that no phase coexistence is allowed here. 

Due to the non-ergodicity of a system, that posesses 2 phases without tunneling in 
between, the trajectory of case (3) is realized in any time dependent adiabatic process 
as used in our correlation dynamics approach for the computation of the effective 
potential. This results in a good convergence of the cluster expansion throughout the 
whole process, because the system always stays in a localized state. For the double-well 
potential with finite potential barrier this is e.g. realized within the ^Q +1 CD(2)-limit, 
i.e. the time-dependent Hartree-Fock method, because at this level of approximation 
tunneling is not yet included. The adiabatic limit of the Qq +1 CD(2) solution then 
simply is the GEP solution, i.e. the lowest energy Gaussian. 

However, for a system with tunneling between the 2 phases we come to the con- 
clusion that the convergence properties of the cluster expansion with respect to the 
adiabatic computation of the lowest energy state are determined by the relative mag- 
nitude of 2 timescales: the specific tunneling time t tun and the time t prop over which 
the system is propagated; the latter quantity in our case is always chosen to ensure suf- 
ficient adiabatic convergence in the single phase region. The tunneling time for a given 
system itself depends on the correlation dynamical truncation scheme; as indicated 
above, t tun is always infinite within the 2-point truncation scheme and the inclusion of 
higher order connected Green functions will in general tend to decrease it. 

If ttun 3> tprop, a dynamical phase mixing will be negligible and our correlation 
dynamics will give useful results for the dynamical time evolution of the system. In 
the two phase region, however, the system will no longer propagate along an adiabatic 
trajectory of lowest energy. Loosely speaking, the correlation dynamical time evolution 
will generate states of lowest energy in the " subspace of sufficiently localized states" , 
the latter of course not being mathematically well defined. 

On the other hand, for t tun <t prop tunneling will lead to a breakdown of correlation 
dynamics in the course of the time integration. This is the reason for the breakdown 
of the $q +1 CT>(2, 3, 4, 5, 6) approach for x / \/JI<3.3 in the upper right part of Fig. |l[ 
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In order to investigate the behaviour of our double-well system with respect to the 
different time scales in a manner directly related to the idealized picture given in Fig. 
D, we now choose a different dynamical process, in which the background field xo is 
changed as a function of time. For the initialization at t = we make use of the 
fact that for xq —>■ oo (strong field limit) the system approaches the classical limit; 
there again the GEP solution can be taken as a starting configuration even though 
the residual interaction is completely taken into account in the time propagation right 
from the beginning. 

We then compute the corresponding external source J as a function of x by impos- 
ing the Hellmann-Feynman theorem (which of course strictly only holds for an exact 
Eigenstate of the system) 

jjE^J) = ±(iB {J)\Hj\MJ)) = (MJ)\^\MJ)) 

= -^ (J)MMJ)) = -MXo(J) , (13) 
where is given by (H). Inserting (g) into ( |T3"D yields 

- /ix (J) = ±{H() = ^j{H 2 ) - /ixo(j) - VJ-^MJ) • (14) 
With 

>■> = (*H &° <j) ) < i5) 

we obtain 

JM = tt-(H 2 ) (16) 
/i ax 

relating the external source J to the effective potential in the usual way. The philosophy 
behind this definition of J(x ) within a dynamical calculation is the interpretation 
already indicated above, that for t tun ^> t prop the system will " quasi-adiabatically" 
evolve in the subspace of localized (single phase) states. 

In Fig. |4]the external source J is plotted as a function of xq/^/I in fg +1 CZ)(2, 3, 4, 5, 6)- 
approximation for a double-well potential with A/(4/i 3 ) = 0.0158, where xo has been 
changed in time according to x (t)/ v //I = x^ tar * / '-JJI—at. The fat solid line is the result 
in $o+iC*i}(2)-approximation with a = 0.00001 c/fm (/j 2 = 1 MeV 2 ), which essentially 
lies on top of the GEP result and is shown for comparison. For large values of a the 
trajectory in the (J — x )-plane can be computed down to x = without any problem 
and exhibits the same qualitative behaviour as the trajectory of case (3) in Fig. |3], 
indicating that for a > 0.01024 c/fm we have t tun 3> t prop ocl/a. 

By reducing a we get into the region where ttun ~ t prop . The inclusion of higher 
order connected Green functions then allows for tunneling, thus leading to a breakdown 
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of the dynamical time propagation due to the dominance of fluctuations in the phase 
coexistence regime. With decreasing a, the value of x , at which this breakdown 
happens, shifts towards the point where J changes sign (xo ~ 3.8). In the exact 
solution this value corresponds to the boundary of the phase coexistence region (cf. 
upper left part of Fig. [1]). 



3 Convergence of the cluster expansion in (1 + 1) 
and (2 + 1) dimensions 

We now turn to the technically much more involved case of (1 + 1)- and (2 + 1)- 
dimensional <3> 4 -theory. Both theories are superrenormalizable, and only a mass renor- 
malization is required. The Lagrangian and the Hamiltonian are given by 



C = -d^d^ - -m 2 $ 2 - 
2 2 4 



(17) 



H 



- I d u x 



1 



n 2 + (V$) 2 + m 2 $ 2 + ^A$ 4 



where m 2 , = m 2 + 5m 2 is the bare mass, m 2 the renormalized mass and v denotes the 
number of spatial dimensions. 5m 2 contains the mass counterterms, which in (2 + 1) 
dimensions are given by the tadpole diagram and the setting sun diagram, and in (1 + 1) 
dimensions by the tadpole diagram only (cf . Jl5|, [L6| ) . 

The field operators and their conjugate momenta are expanded into plane waves in 
a box with periodic boundary conditions according to 



= E 



u(x) = £ n a ^o 



x 



(19) 



with ip a (x) = -j^e lkaX ; V = L u and k a chosen according to the boundary conditions. 

For our numerical simulations in (1 + 1) dimensions we use L = 100 fm, m = 
10 MeV and the 21 lowest lying plane waves; in (2 + 1) dimensions we use L = 20 fm, 
m = 10 MeV and the 29 lowest lying states. For further details of the explicit set of 
equations solved we refer the reader to refs. |15|, [16 . 

In Fig. H we show the quantities 



E<*«*»>c/E<*«**> 



a/3 



a/3 



<$ 2 (*)> c /($ 2 (*)> 



in 



2) 



(20) 



E<*«*/>*7>e/E<*«*/»*7> 

a/3-y a/37 



($ 3 (x)) c /($ 3 (x)) 



in 



3) 



(21) 
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and 



E 

af3ji 



a/375 



(<& 4 (z)) c /<$ 4 (a;)> 



in 



4) 



(22) 



as a function of $0 — (&( x )) — L u / 2 (Q(k = 0)). Note that due to translation invari- 
ance all Green functions with J2 ai k ai ^ have to vanish. The quantities (20"D, ( pTj ) 
and (^) have been chosen as a measure for the relative importance of the higher order 
connected Green functions in analogy to (^|). All results have been obtained within the 
$ 4 +1 CD(2, 3, 4)-approximation (i.e. including everything up to the connected 4-point 
function), which is the lowest order in correlation dynamics to which the setting sun 
diagram contributes, which diverges logarithmically in (2+1) dimensions (cf. flEfl). 
The coupling has been switched on in time according to 



A 



4m 3 ~ 



fit 



(23) 



with the free (perturbative vacuum) solution as initialization. This method is more 
efficient than switching on the residual interaction with the GEP solution as an initial- 
ization, since it yields results for the lowest energy state over a whole range of coupling 
constants in only a single time dependent run. In the (0 + l)-dimensional case, the 
method of switching on the residual interaction (leading to a correlated state for only 
one value of the coupling constant) had to be chosen, because the double well potential 
has no "unperturbed phonon state", because for zero coupling the potential is a har- 
monic oscillator with negative mass and thus not bounded from below. The vacuum 
expectation value $0 has been split off as a classical background field and the 1-point 
function has not been propagated explicitly in analogy to the treatment of x in the 
(0 + l)-dimensional case. The value j3 = 0.05 c/fm has been chosen small enough in 
order to observe an adiabatic behaviour of the time dependent solution. 

The results for (1 + 1) dimensions are shown in the left column and those for (2 + 1) 
dimensions in the right column of Fig. |5] for various values of the coupling constant. 
As in Fig. [I], the peaks result from taking the absolute value of a quantity that changes 
sign and have no physical meaning. In both, (1 + 1) and (2 + 1) dimensions, the relative 
importance of the connected 3- and 4-point functions increases with increasing coupling 
and decreasing $0; however, for all parameters considered this relative importance 
decreases in going from n to n + 1, and moreover, the highest contribution of the 
connected 3-point function to the full 3-point function is only about 20 %. The cluster 
expansion thus shows satisfactory convergence throughout the investigated parameter 
range. 

The critical couplings for a transition into the symmetry broken phase in (1 + 1) and 
(2 + 1) dimensions within the $ 4 +1 C-D(2, 3, 4)-approximation are (for the present spe- 
cific choice of numerical parameters specifying the plane wave basis set) A /4m 2 = 2.247 
and A/4m = 0.369, respectively |TR TB\. Below these couplings a good convergence of 
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the cluster expansion can be expected, since the system is in a single phase, i.e. it can 
be described by a localized wave functional. The fact, that we still observe this con- 
vergence above the critical coupling indicates that the time t tun required for building 
up a two phase configuration via tunneling within the $£ +1 CD(2, 3, 4)-approximation 
is infinite or at least large compared to the time t prop needed for the " quasiadiabatic" 
propagation of the system along the lowest energy localized state trajectory in the sense 
of Sect. H In the latter case, the system would (over a certain range of t prop oc 1//3) 
seem to converge against an asymptotic configuration as t prop increases, and then sud- 
denly change its behaviour when tunneling kicks in. We have only been able to observe 
such a behaviour in the (0 + 1) dimensional case, but within the investigated range of 
f3- values not in (1 + 1) and (2 + 1) dimensions. 

In fact, it is well known that tunneling is suppressed as the number of space-time 
dimensions increases. In (1 + 1) dimensions, by means of a soliton-antisoliton pair, it 
can be shown that tunneling between the two vacua with $o = ±$™ m is still possible 
J?|. In (2 + 1) dimensions, however, no tunneling is expected [[7|]; this would be the 
lowest space-time dimension with real ground state symmetry breaking, i.e. with an 
effective potential which is flat between — $™ n and +$' nm . 

A simple qualitative argument [[?[] pointing towards the suppression of tunneling in 
higher space-time dimensions is the following: Consider a bubble of phase 1 within a 
region of phase 2. The surface region of this bubble then has a higher energy density 
than the surrounding. In (1 + 1) dimensions, the surface region is given by the envi- 
ronment of the two end-points of the bubble interval (the soliton and the antisoliton) 
and its size is therefore independent of the size of the bubble; consequently the bubble 
can expand without increasing the total energy of the system. Generally, in (d + 1) 
dimensions, the size of the surface region behaves as so that for d > 2 the bub- 

ble cannot expand without a corresponding increase in the total energy density. The 
above geometrical considerations are similar to those used for the proof of Derrick's 
theorem |23| (non-existence of topological solitons in scalar theories with d > 2), where 
a rescaling of the classical field according to $(x) — > $(Ax) is used to show that any 
potential topological soliton solution would be energetically unstable against variation 
of A. 

4 Summary 

In this work we have explored the convergence properties of the cluster expansion 
approach and presented criteria within the theory itself, that allow to conclude about 
its convergence also in those cases, where the exact solution is not known. As an 
example we have studied scalar theories with quartic self coupling in (0 + 1), (1 + 1), and 
(2 + 1) dimensions. We have emphasized the aspect of a nonperturbative computation 
of the lowest energy state with given field expectation value, since our previous work 
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15 , [16] aimed at the evaluation of the effective potential. 

The $ 4 theory in (0 + 1) dimensions is equivalent to a quantum mechanical an- 
harmonic oscillator for positive squared mass and a quantum mechanical double-well 
potential for negative squared mass. In the case of an anharmonic oscillator with 
A/(4m 3 ) = 10 we obtained a good convergence of the cluster expansion as well as an 
excellent agreement between the exact solution and the correlation dynamical solution; 
this could be explained by the fact that the system is in a single phase (localized state) 
configuration. 

In case of an extreme double-well potential with A/(4/i 3 ) = 0.0158 the cluster 
expansion in the exact solution breaks down as soon as the ground state for a given 
expectation value (x) can approximately be described by a Maxwell construction, i.e. 
by a coherent superposition of two single phase configurations. The shrinking of the 
region with dominant fluctuations within the correlation dynamical calculation could 
be traced back to the fact, that the time t prop , chosen according to the requirement of 
an adiabatic propagation in the single phase region, is small compared to the tunneling 
time t tun of the system for certain values of (x). 

Our analysis for the QFT systems in (1 + 1) and (2 + 1) dimensions demonstrated, 
that even in the phase coexistence (non-localized state) region, i.e. for coupling con- 
stants higher than the critical coupling for ground state symmetry breaking, the cor- 
relation dynamical approach yields a good convergence of the cluster expansion. Here, 
the tunneling times within our dynamical calculation are infinite or at least very large 
compared to the time required for an adiabatic propagation; this behaviour is traced 
back to the fact that tunneling is suppressed in (1 + 1) dimensions and probably absent 
in (2 + 1) dimensions. 

The main result of this article can therefore be summarized as follows: the appli- 
cability of correlation dynamics (which is based on an expansion in terms of connected 
Green functions) to the adiabatic computation of the lowest energy state with given 
field expectation value is given if a) the exact stationary solution is a single phase, 
and b) if the required propagation time for adiabatic or " quasiadiabatic" (cf. Sect. ^|) 
convergence is small compared to the tunneling time. This is in fact the case in higher 
space-time dimensions since there tunneling is progressively suppressed. This is not 
only important for lowest energy state calculations as carried out in this article, but 
also for applications of correlation dynamics to the nonequilibrium dynamics of systems 
with phase transitions. Furthermore, the relative ratios (|20| - |2"2"|) provide an intrinsic test 
of the correlation dynamical method also for those field theoretical problems, where 
the exact solution is not known. 
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A The cluster expansion 



The explicit form of the cluster expansion can be derived from the generating func- 
tional of full and connected Green functions (cf. e.g. [0), Z[J,a] and W[J, a], given 
by 

Z[J, a] = Tr {pT [e l I ^ + M^)^)+-(z)n(£))] j ^ 
and 

Z[J,a] = e w[J > a] , (25) 

respectively, where T is the time ordering operator, p is the statistical density operator 
describing the state of the system (Trp = 1) and x denotes (x,t). We start with the 
cluster expansions for the time-ordered Green functions: 

(*(x)> = (Q(x)) e , (II(*)> = (U(x)) c , (26) 

(T$(x 1 )$(.x 2 )) = Km -^—-^—e w ^ 

= lim -7^l(-7Tr^ W l J ^]) eW[J ' a] ) 
j^oiSJ( Xi ) [\i5J(x 2 ) J J 

= lim ( f . * r— ^ r W[J, a]) + ( rJr^W[J, a]) ( rjJ-^W[J, a}) ) e w ^ 
J,<t->o{\i8J(xi)i8J(x 2 ) 7 Vz5J(xi) 7 W^s) 7 J 

= (T$(x 1 )$(x 2 )) c + ($(x 1 ))($(x 2 )) , (27) 
where (-) c denotes the connected part of the expectation value. Analogously we obtain 

(Tn(zi)$(x 2 )) = (Tn(xx)$(x 2 )) c + (n(xi))($(x 2 )) , 
(Tn(z!)n(z 2 )) = (rn(x!)n(x 2 )) c + (n(xx))(n(x 2 )) , 

+(T$(x 1 )$(x3)) c ($(x 2 )} + (r$(x 2 )$(x3)) c ($(x 1 )} + ($(x 1 )}($(x 2 ))($(x 3 )) , ... . 

(28) 

The expressions for equal-time Green functions are obtained by taking the well-defined 
equal-time limit which yields the appropriate operator ordering in the cluster expan- 
sions. We arrive at 

(*(*)) = ($(*)) e , (n(*)> = (n(x)) c , 
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($(ar 1 )$(ar 2 )) = ($(x 1 )$(x 2 )) c + <$(a:i))<$(z 2 )) , 

(n(x!)$(x 2 )> = (n(x 1 )$(x 2 )) c + (n(x!))($(x 2 )) , 
(n(xOn(x 2 )> = (n^n^)}^^^)}^)) , 

($(xi)$(x 2 )$(x 3 )) = ($(xi)$(x 2 )$(x 3 )) c + ($(xi)$(x 2 )) c ($(x 3 )) 

+ ($(X 1 )$(X3)} C ($(X 2 )} + ($(X 2 )$(X3)} C ($(X 1 )} + ($(X 1 )}($(X 2 )}($(X 3 )) , - , 

(29) 

where all (equal) time arguments have been suppressed. In view of their length the 
cluster expansions for the other Green functions required for our calculations are not 
explicitly given here, but can be found in The explicit forms of the cluster expan- 
sions of (x"), n — 1, • • • , 6 (required for our analysis in (0 + 1) dimensions) are given 
in appendix |B|. 



B Cluster expansion of (x n ) 

Applying the derivation of appendix we obtain the following expressions for the 
cluster expansion of (x n ), n — 1, • • • , 6 in $Q +1 -theory using $(x) — > x: 

(x) = (x) c , (30) 

<x 2 > = <*>c + <x 2 > c , (31) 

(x 3 ) = (x)^ + 3(x) c (x 2 ) c +(x 3 ) c , (32) 

<x 4 > = <x> 4 + 6(x) 2 (x 2 ) c + 3(x 2 ) 2 + 4(x) c (x 3 ) c + (x 4 ) c , (33) 

(x 5 ) = (x)« + 10(x) 3 (x 2 ) c + 15(x) c (x 2 ) 2 + 10(x) 2 (x 3 ) c + 10(x 2 ) c (x 3 ) c 

+5(x) c (x 4 ) c +(x 5 ) c , (34) 

(A = (*}t + 15(x) 4 (x 2 ) c + 45(x) 2 (x 2 ) 2 + 15(x 2 ) 3 + 20(x) 3 (x 3 ) c 
+60(x) c (x 2 ) c (x 3 ) c + 10(x 3 ) 2 + 15(x) 2 (x 4 ) c + 15(x 2 ) c (x 4 ) c + 6(x) c (x 5 ) c + (x 6 ) c . (35) 
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C The GEP approximation in (0 + 1) dimensions 



We consider the Hamiltonian 

H = l - V 2 + ^m 2 x 2 + ^Ax 4 = X - V 2 + X -m 2 {x + x ) 2 + -\{x + x ) 4 , (36) 

where m 2 can be positive or negative. The GEP (Gaussian effective potential) in (0 + 1) 
dimensions is given by a constrained variational calculation with a Gaussian trial wave 
function, i.e. by the ansatz 




This ansatz can then be seen as the ground state of a variational harmonic Hamiltonian 
of the form 

H= \ p2 + (38) 

with an effective mass M. The connection to the field theoretical Hartree-Fock- 
Bogoliubov ansatz becomes obvious by noting, that Q57| ) is the quasiparticle vacuum 
with respect to Bogoliubov-transformed creation and annihilation operators. Keeping 
this last point in mind, it is also straightforward to see why in a Gaussian trial state 
all connected Green functions of higher order than the two-point function vanish; this 
property follows directly from Wick's theorem. 
The resulting equation for the effective mass is 

M 2 - m 2 = 3A (x 2 + JL) , (39) 

which is equivalent to the Gap equation following from the Bogoliubov- ansatz. 
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Figure Captions 



fig.l: The ratios |(x n ) c /(x n )| for n=2,...,6 as a function of xq/^JI = (x)/y / /7 for 
the double-well potential with A/(4/i 3 ) = 0.0158 (upper two pictures) and 
for the anharmonic oscillator potential with A/(4m 3 ) = 10 (lower two pic- 
tures). L.h.s.: exact results; r.h.s.: results within the $g +1 CD(2, 3, 4, 5, 6)- 
approximation obtained by time dependently switching on the residual in- 
teraction. 

fig. 2 : Absolute squares of the ground state wave function ^o(x) and the wavefunc- 
tion of the first excited state ^i(x) (both obtained from the exact solution) 
as a function of x and of f(J) = (\& | x |^o) i n the case of a double-well 
potential with A/(4/z 3 ) = 0.0158 and an external source J. 

fig. 3: Qualitative picture of the external source J as a function of x = (x) for 
3 different cases: (1) no symmetry breaking, (2) Maxwell construction in 
case of symmetry breaking, (3) analytical continuation into the two phase 
region of the Maxwell contruction in case of symmetry breaking. 

fig. 4: External source J as a function of xq/^/JI in Qq +1 CD(2, 3, 4, 5, 6)-approximation 
for a double-well potential with A/(4/x 3 ) = 0.0158; xo is changed in time 
according to x (i)/ v //I = x^ art / \fji — at. The fat solid line is the result in 
$^ +1 CL'(2)-approximation with a = 0.00001 c/fm (/i 2 = 1 MeV 2 ). 

fig-5 : The ratios \j2 a p(®a®p)c / E a p(®«®p) \ (n=2), |E a ^(^a^ 7 )c / £a/3 7 ($a$/3$ 7 ) 

{n = 3) and |£ a/ ^($ a $/j$ 7 $«)c / £ a/ ^($ a $/9$ 7 $«) | ( n = 4) as a func- 
tion of $o — (*&( x )) f° r (1 + 1) (1-h.s.) and (2 + 1) (r.h.s.) dimensions in 
<£> 4 C.D(2, 3, 4)-approximation; all results have been obtained by switching 
on the coupling adiabatically. 
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